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HMXB, ULX and star formation 
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Based on recent X-ray observations of the Milky Way, Magellanic Clouds and nearby starburst galaxies we study 
population of high mass X-ray binaries, their connection with ultra- luminous X-ray sources and relation to the star 
formation. Although more subtle SFR dependent effects are likely to exist, the data in the \g{Lx) ~ 36 — 40.5 
luminosity range are broadly consistent with existence of a "universal" luminosity function of HMXBs, which 
can be roughly described as a power law with differential slope of ~ 1.6 and a cutoff at lg{Lx) ~ 40.5. The 
ULX sources found in many starburst galaxies occupy the high luminosity end of this single slope power law 
distribution, whereas its low luminosity part is composed of "ordinary" high mass X-ray binaries, observed, e.g. 
in the Milky Way and Magellanic Clouds. 

As the normalization of the "universal" luminosity function is proportional to the star formation rate, the 
number and/or the collective X-ray luminosity of HMXBs can be used to measure the current value of SFR in 
the host galaxy. Distant (unresolved) starburst galaxies observed by Chandra at redshifts of « ~ 0.2 — 1.3 obey 
the same Lx—SFR relation as local galaxies, indicating that the ULXs at these redshifts were not significantly 
more luminous than those found in nearby galaxies. 



1. INTRODUCTION 

An unusual class of compact sources - ultra- 
luminous X-ray sources, has been discovered in 
nearby galaxies more than a decade ago [2,3]. 
Although bright, Lx > 10"^^ erg/s, point-like 
sources are found both in young star forming 
galaxies and in old stellar population of elliptical 
galaxies, the most luminous and exotic objects 
are associated with actively starforming galaxies. 
Their nature and relation to more ordinary X-ray 
binaries is still a matter of a significant debate. 
Based on a simple Eddington luminosity argu- 
ment, they appear to be powered by accretion 
onto an intermediate mass object - a black hole 
with the mass in the hundreds-thousands solar 
masses range. However, a number of alternative 
models have been considered as well - from coUi- 
mated radiation to ^stellar mass black holes, rep- 
resenting the high mass tail of the standard stellar 
evolution sequence and accreting in the near- or 
slightly super-Eddington regime. 

Sub-arcsec angular resolution of Chandra ob- 
servatory opened a new era in studying popula- 
tion of compact sources in nearby galaxies. For 



the first time an opportunity was presented to 
observe compact sources in a nearly confusion 
free regime and to investigate their relation to 
fundamental parameters of the host galaxy, such 
as its stellar mass and star formation rate. De- 
pending on the mass of the optical companion. 
X-ray binaries are subdivided in to two classes - 
low and high mass X-ray binaries, having signifi- 
cantly different evolutionary time scale, ^ 10^^^ 
and ^ 10^~^° years respectively [12]. The prompt 
emission from HMXBs makes them a potentially 
good tracer of the very recent star formation ac- 
tivity in the host galaxy [11]. The LMXBs, on the 
other hand, have no relation to the present star 
formation, but, rather, are related to the stellar 
content of the host galaxy [4]. Chandra obser- 
vations of the nearby galaxies present a possibil- 
ity to verify this simple picture and to calibrate 
the relation between the HMXB population and 
the current value of SFR. This opens a new way 
to determine the star formation rate in relatively 
closeby galaxies as well as in young and very dis- 
tant ones. Detailed study of the population of the 
compact sources in star forming galaxies can also 
shed the light on the mistcry of the ULXs. 
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Figure 1. Left: The luminosity functions of compact X-ray sources in nearby galaxies obtained by 
Chandra. Right: The luminosity functions for the same galaxies scaled by the ratio of their star formation 
rate to that of Antennae. 



2. "UNIVERSAL" 



LUMINOSITY 



FUNCTION OF HMXBs 

In order to study the HMXB population and its 
relation to star formation we selected a number of 
nearby late type galaxies observed by Chandra, 
based primarily on two criteria, (i) The galaxy 
can be spatially resolved by Chandra, so that the 
contribution of a central AGN can be discrimi- 
nated and the luminosity functions of the com- 
pact sources can be constructed without severe 
confusion effects, (ii) The galaxy has sufficiently 
high SFR/mass ratio, to ensure that the popula- 
tion of X-ray binaries is dominated by HMXBs 
and the LMXB contribution of can be safely ig- 
nored [5]. For the Milky Way we explicitly se- 
lected HMXBs, based on results of [7]. For each 
galaxy, the star formation rate was determined 
combining the results from conventional SFR in- 
dicators, [9] (FIR, UV, H„ and radio). The SFR 
values in the sample range from ~ 0.15 to ~ 7 
Mq/ji. Further details and references to the X- 
ray data are given in [8]. 

The Fig.l (left panel) shows observed luminos- 



ity functions. They are characterized by a large 
spread in the number of sources and in the lu- 
minosity of the brightest source. However, after 
being rescaled to the same value of the star forma- 
tion rate, they appear to match each other both 
in the slope of the distribution and the normal- 
ization (right panel in Fig.l). Although a finite 
dispersion might still exist, especially at the high 
luminosity end, the rescaled luminsity functions 
occupy a rather narrow band in the N — Lx plane, 
despite of large dynamical range of the star for- 
mation rates, a factor of ~ 50. We further com- 
bine the data for five nearby starburst galaxies 
with the best known luminosity functions, having 
cumulative star fromation rate of « 16 Mq/yt, 
and compare them with the luminosity distribu- 
tions of two low SFR galaxies, the Milky Way and 
SMC (Fig. 2). As previously, a good agreement 
both in the slope and normalization of luminos- 
ity distributions is obvious. The fit to combined 
data with a power law distribution gives the slope 
of a w 1.6, a cut-off at icut « (2 - 3) • lO""^ erg/s 
and the normalization proportional to the star 
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Figure 2. Left: Combined luminosity function of compact X-ray sources in star burst galaxies and the 
luminosity functions of HMXBs in the Milky Way and SMC. The thick grey line is the best fit to 
the starburst galaxies, Eq.(l). Right: Average differential luminosity function of the galaxies from our 
Chandra sample. The straight line is a power law defined by eq.(l). The shaded area illustrates the 
amplitude of systematic errors (90% confidence level) associated with uncertainties in the adopted SFR 
values (assuming 30% relative error) and the distances (20% relative error). 



formation rate: 

^ = (3.3111) ■ SFR- L^i-^^^^-^' (1) 

where SFR is formation rate for massive stars, 
M > 5Mq. Within the accuracy of the present 
analysis, this can be regarded as a "universal" lu- 
minosity function of high mass X-ray binaries in 
nearby star forming galaxies. The Fig. 2 shows 
the average differential luminosity function of 
HMXBs obtained using the data of all galaxies 
from our Chandra sample. 

3. HIGH MASS X-RAY BINARIES AS A 
STAR FORMATION INDICATOR 

As the population of high mass X-ray binaries 
(i.e. the normalization of the luminosity func- 
tion and the number of sources) is proportional 
to SFR, eq.(l), they can be exploited to measure 
star formation rate in the host galaxy. Such a 



method, based on the X-ray emission of a galaxy, 
might circumvent one the main sources of uncer- 
tainty of conventional SFR indicators - absorp- 
tion by dust and gas [9]. Indeed, galaxies are 
mostly transparent to X-rays above ~ 2 keV, ex- 
cept for the densest parts of the most massive 
molecular clouds. 

For nearby, spatially resolved galaxies, the star 
formation rate can be determined via a direct fit- 
ting of the luminosity distribution of the compact 
sources or counting the number of sources above 
a given luminosity threshold. For distant, unre- 
solved galaxies one can use the total luminosity 
of the galaxy. 

There are, however, two main complicating fac- 
tors in using the X-ray luminosity of a galaxy as a 
SFR indicator: (i) X-ray emission from the cen- 
tral supermassive black hole, which even in the 
case of a low luminosity AGN can easily out- 
shine X-ray binaries. This does not present a 
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Figure 3. Left: The probability distributions p(Ltot/SFR) for different values of SFR. The vertical dashed 
line shows the expectation mean, defined by eq.(2). Right: The ix— SFR relation. The open circles are 
nearby galaxies observed by Chandra, the filled triangles are nearby galaxy observed by ASCA and 
BeppoSAX, for which only total luminosity is available, the filled circles are distant star forming galaxies 
from the Hubble Deep Field North. The thick grey line is relation between the star formation rate and 
the most probable value of the total luminosity, predicted from the "universal" luminosity function of 
HMXBs, the shaded area shows 67% intrinsic spread of the Lx— SFR relation, the dashed line is the 
expectation mean, defined by eq.(2). 



problem in nearby {D < 20 — 30 Mpc) galaxies 
as the Chandra resolution is sufficient to sepa- 
rate contribution of the central source, but might 
become an important source of contamination in 
spatially unresolved galaxies, (ii) Contribution 
of low mass X-ray binaries. This equally affects 
nearby and distant galaxies, as LMXBs can not 
be easily separated from HMXBs using the X- 
ray data in the Chandra bandpass, and the op- 
tical identifications are (potentially) available for 
the most nearby galaxies only. However, as the 
LMXB population scales linearly with the stellar 
mass of the host galaxy [4] , their contribution can 
be estimated and, potentially, corrected for. Ob- 
viously, at sufficiently high values of SFR/M, it 
can be safely neglected, [5] . 

To further calibrate the Lx— SFR relation and 
verify its applicability to distant galaxies, we ex- 



tended our sample to include the spatially un- 
resolved galaxies observed by ASCA and Bep- 
poSAX and distant star forming galaxies detected 
by Chandra in the Hubble Deep Field North [1] 
(see [8] for details and references). As for the 
galaxies from our spatially resolved Chandra sam- 
ple, star formation rates were determined by the 
conventional, non X-ray methods. The combined 
data are shown in the L^— SFR plane in Fig. 3. 

3.1. Effects of statistics and ix — SFR re- 
lation 

A seemingly obvious expression for the total 
luminosity can be obtained integrating the lumi- 
nosity distribution (1): 

(itot) = J —LdLo^SFR (2) 
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implying, that the total luminosity is propor- 
tional to the star formation rate. We note, how- 
ever, that the quantity of interest is a sum of the 
luminosities of discrete sources: 

Ltot = ^ Lk (3) 

fc 

with Lfc obeying a power law probability distribu- 
tion given by eq.(l). For a sufficiently flat slope, 
a < 2, the total luminosity Ltot will be defined 
by the brightest sources, corresponding to the 
high luminosity end of the power law distribution 
eq.(l). Their actual number in a galaxy will obey 
the Poisson distribution P(n,/Lt) = /Lt"e^^/n! 
(with n ocSFR), which for small values of /i is 
significantly asymmetric. As a consequence, for a 
galaxy with small SFR, the probability distribu- 
tion p(Ltot) will be also strongly asymmetric, as 
illustrated by the left panel in Fig. 3. Because of 
its skewness, the mode of the p(Ltot) distribution 
- the value of Ltot that would be most likely mea- 
sured in an arbitrarily chosen galaxy, is not equal 
to the expectation mean defined by eq.(2). Only 
in the large SFR limit, when there are sufficiently 
many sources with luminosities L ~ Lcut, these 
two quantities become close to each other. ^ 

The difFcrcnce between these two quantities is 
further illustrated in the right panel in Fig.3, 
showing the predicted Lx-SFR relation, calcu- 
lated using the parameters of the "universal" 
HMXB luminosity function, eq.(l). The solid line 
in the figure shows the SFR-dependence of the 
mode of the probability distribution p(Ltot) and 
predicts the most probable value of the X-ray lu- 
minosity of a randomly chosen galaxy. If observa- 
tions of many (different) galaxies with close values 
of SFR are performed, the obtained values of Ltot 
will obey the probability distribution depicted in 
the left panel of Fig.3. The average of the mea- 
sured values of Ltot will be equal to the expec- 
tation mean given by eq.(2) and shown by the 
dashed straight lines in the left and right panels 
of Fig.3. Due to the properties of the probability 
distribution p(Ltot) these two quantities are not 
identical in the low SFR limit, when the total lu- 

^ Obviously in the case of e.g. flat (dN/dL =const) or 
Gaussian flux distribution the most probable value of Ltot 
always equals to the expectation mean defined by eq.(2). 



minosity is defined by a small number of the most 
luminous sources. 

3.2. Lx— SFR relation: predicted vs. ob- 
served 

The right panel in Fig.3 compares the data with 
the predicted L^— SFR relation. Good agree- 
ment, both in the non-linear low SFR regime and 
at high SFR values is apparent. The relations be- 
tween total X-ray luminosity of a galaxy due to 
HMXBs and the star formation rate are: 

in the linear regime and 

in the non-linear regime, corresponding to SFR< 

4.1 Mq yr-i (L2-10 kcv < 2.8 • 10^" erg s'l) 
The former, linear, regime of the Lx~SFR re- 
lation was studied independently by [10] based 
on ASCA and BcppoSAX data. Note that their 
equation (12) agrees with our eq.(4) within ~ 30 
per cent. 

Due to skewness of the probability distribution 
p(Ltot), large and asymmetric dispersion around 
the solid curve in Fig.3 is expected in the non- 
linear low SFR regime. This asymmetry is al- 
ready seen from the distribution of the points in 
Fig.3 - at low SFR values there are more points 
above the solid curve, than below. Moreover, the 
galaxies lying significantly above the solid and 
dashed curves in Fig.3 should be expected at low 
SFR and will inevitably appear as the plot is pop- 
ulated with more objects. Such behavior difi'ers 
from a typical astrophysical situation and should 
not be ignored when analyzing and fitting the 
Ljsf-SFR relation in the low SFR regime. In par- 
ticular, the standard data analysis techniques - 
least square and fitting become inadequate. 

4. DISCUSSION 

4.1. High luminosity cut-off in the HMXB 

luminosity function 
The existence of the linear regime in the Lx- 
SFR relation is a direct consequence of the cut- 
off in the luminosity function. The position of 
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Figure 4. Dependence of the SFR relation 
on the maximum luminosity of the sources, Lcut- 
The three curves, corresponding to different val- 
ues of Lcut coincide in the non-linear low SFR 
regime but differ in the position of the break be- 
tween linear and non-linear regimes. The data 
are the same as in Fig. 3. 



the break between non-linear and linear parts of 
the Lx^SFR relation depends on the slope of the 
luminosity function and the value of the cut-off 
luminosity (Fig.4): SFRbroak oc L'^~t^. This al- 
lows one to constrain parameters of the luminos- 
ity distribution of compact sources using the data 
of spatially unresolved galaxies. 

Agreement of the predicted SFR relation 
with the data both in high and low SFR regimes 
gives an independent confirmation of the exis- 
tence of a cut-off in the luminosity function of 
HMXBs at Lcut ~ several x 10'*° erg/s (Fig.3,4). 
It also confirms that Lx— SFR data, including 
the high redshift galaxies from Hubble Deep Field 
North, are consistent with the HMXB luminosity 
function parameters, derived from significantly 
fewer galaxies, than plotted in Figs. 3, 4. 



4.2. Ultraluminous X-ray sources 

One of the surprising results of this study is 
a smooth, single slope power law shape of the 
average luminosity function of compact sources 
in star forming galaxies, without any significant 
steps and features in a broad luminosity range, 
lg{Lx) ~ 36 - 40.5 (Fig.2). The high luminos- 
ity end, Ig(Lx) > 39, of this distribution corre- 
sponds to ultraluminous X-ray sources. Its low 
luminosity end, on the other hand, is composed 
of ordinary X-ray binaries, powered by accretion 
onto a ~stellar mass compact objects. This re- 
sult constrains the range of possible models for 
ULXs. Their frequency and luminosity distri- 
butions should be a smooth extension towards 
higher luminosities of that of "ordinary" ~stellar 
mass systems, emerging from the standard stellar 
evolution sequence. Although some of the ULXs 
might be indeed rare and exotic objects, it ap- 
pears that majority of them cannot be a com- 
pletely different type of the source population, 
but, rather, represent the high mass, high M tail 
of the HMXB population. 

The luminosity of ULXs in the nearby galaxies 
has a maximum value of the order of \g{Lx) ~ 
40.5. The fact that the galaxies from the Hub- 
ble Deep Field North obey the same L^— SFR 
relation (Fig.4), implies, that the ULXs at the 
redshift of z '--^ 0.2 — 1.3 were not significantly 
more luminous, that those observed in the nearby 
galaxies. 

4.3. Intermediate mass black holes 

The hypothetical intermediate mass black 
holes, probably reaching masses of ^ 1O^~^M0, 
might be produced, e.g. via black hole merges 
in the dense stellar clusters, and can be associ- 
ated with extremely high star formation rates. 
To accrete efficiently, they should form close bi- 
nary systems with normal stars or be located in 
dense molecular clouds. It is natural to expect, 
that such objects are significantly less frequent 
than '--^stellar mass black holes. The transition 
from the ~stellar mass BH HMXB to interme- 
diate mass BHs should manifest itself as a step 
in the luminosity distribution of compact sources 
(Fig. 5, left panel). If the cut-off in the HMXB 
luminosity function, observed at Ig(Lcut) ^ 40.5 
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Figure 5. Illustration of the effect of hypothetical intermediate mass black holes on the Ljf— SFR rela- 
tion. Left: The luminosity function of compact sources at different levels of star formation rate. Right: 
Corresponding Lj^— SFR relation. The thin straight line shows the linear relation. 



corresponds to the maximum possible luminosity 
of ~stellar mass black holes and if at L > Lcut 
a population of hypothetical intermediate mass 
BHs emerges, it should lead to a drastic change 
in the slope of the Lx^SFR relation at extreme 
values of SFR [6] (Fig. 5, right panel). Therefore, 
observations of distant star forming galaxies with 
very high SFR might be an easy way to probe the 
population of intermediate mass black holes. 
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